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Abstract 

The  future  of  the  Army’s  air  delivery  mission  includes 
the  use  of  precision-guided  autonomous  airdrop 
methods  to  resupply  troops  in  the  field.  High-glide 
systems,  ram-air  parafoil-based,  allow  for  a  safe 
standoff  delivery  as  well  as  wind  penetration.  This 
paper  addresses  the  development  of  a  six-degree-of- 
freedom  model  of  a  low-aspect  ratio  controllable 
parafoil-based  delivery  system.  The  model  is  equally 
suitable  for  modeling  and  simulation  and  for  the  design 
of  guidance,  navigation  and  control  (GNC)  algorithms. 
This  gliding  parafoil  model  was  developed  in  the 
MATLAB/Simulink®  environment.  Apparent  mass  and 
inertia  effects  are  included  in  the  model.  Initial  test 
cases  have  been  run  to  check  model  fidelity. 


Introduction 

The  tempo  of  the  battlefield  in  today’s  Army  is 
extremely  fast-paced.  The  Desert  Storm  war 
demonstrated  the  rapid  pace  at  which  a  land  battle  can 
be  executed.  The  Army’s  transformation  initiatives  are 
geared  towards  increasing  the  ability  of  the  Army  to 
respond  rapidly  and  improve  the  agility  and  mobility  of 
its  forces.  A  key  part  of  this  transformation  is  the  ability 
to  get  materiel  where  it  needs  to  be,  when  it  needs  to  be 
there.  In  addition,  the  US  Air  Force  desires  to  improve 
the  survivability  of  its  air  delivery  aircraft  by  increasing 
the  ability  to  drop  payloads  from  higher  altitudes  than 
currently  dropped  from  and  further  from  the  drop 
zones.  Therefore,  the  services  require  the  ability  to 
deliver  payloads  from  high  altitude,  with  large  standoff 
to  achieve  precision  accuracies  from  the  desired  impact 
point. 
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Numerous  papers  in  the  area  of  parafoil  system  design 
have  been  published  over  the  past  50  years.  They  cover 
the  following  major  topics: 

■  Advanced  computational  methods  for  the  3D 
flow  simulation  around  the  parafoil  canopy;1"3 

■  Wind-tunnel  experiments;4"7 

■  Real  drop  experiments.8"18 

Although  each  of  these  papers  topics  addresses  the 
problem  differently,  they  all  share  an  important  feature 
-  verification  of  corresponding  mathematical  models 
using  real  drop  data. 

Despite  the  50-year  research  effort,  the  aerodynamics  of 
a  flexible  body  such  as  a  parafoil  canopy  still  contains 
some  unknowns  and  uncertainties.  For  instance, 
because  of  the  flexible  nature  of  a  parafoil's  canopy, 
even  aerodynamic  experiments  in  wind  tunnels4  may 
provide  variable  data  for  both  longitudinal  and  lateral 
channels,  not  explicitly  applicable  for  the  modeling. 

Analysis  of  previous  work  ’  ’  that  addresses 
model  verification  leads  to  the  conclusion  that  major 
differences  between  them  lie  in  the  way  authors  account 
for  the  influence  of  numerous  interrelated  parameters 
on  the  motion  of  entire  system.  This  paper  does  not 
intend  to  discuss  various  numerical  techniques  of 
parameter  identification  (PID),23  but  briefly  mentions 
the  physical  issues  (multicriteria  essence)  behind  the 
identification  process.  This  paper  employs  the  same 
methodology  authors  have  already  used  for  the  PID  of  a 
controllable  circular  parachute.24,25  The  key  feature  of 
this  original  technique  is  to  separate  the  influence  of 
different  dominant  factors  (apparent  mass  terms, 
aerodynamic  coefficients)  by  employing  different 
portions  of  the  flight  test  data. 

The  paper  is  organized  as  follows.  The  background 
section  of  the  paper  reviews  existing  parafoil  models 
and  discusses  several  modeling  issues  not  completely 
resolved  to  date.  Formulation  of  the  problem  follows 
next.  The  paper  then  proceeds  with  a  development  of 
the  mathematical  model  of  a  double-skin  parafoil.  This 
includes  the  basic  equations  of  motion  in  matrix  form 
convenient  for  further  implementation  within  the 
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M  athworks  M  ATLAB/Simulink  environment.  The 
apparent  mass  model  is  then  presented.  The  following 
section  provides  analysis  and  computation  of  the 
aerodynamic  forces  and  moments.  The  paper  then 
contains  examples  of  simulations  including  system 
response  on  certain  control  patterns.  Preliminary  results 
of  the  model's  parameters  identification  based  on 
comparison  with  the  airdrop  data  available  are  also 
introduced.  The  paper  ends  with  the  summary  of 
obtained  results. 


C  onfiquration  of  the  Descending  System 

The  autonomously  guided  Pegasus  delivery  system  is 
comprised  of  a  ram-air  canopy,  a  GPS-based  guidance 
and  control  unit,  and  payload.  A  photo  of  the  system  is 
shown  in  Figure  1. 


Fig.  1  Pegasus  parafoil  (from  Ref.  26) 


The  main  canopy  is  an  eleven-cell  ram-air  design.  It  is 
made  of  1.1-ounce  F-lll  nylon  with  cotton 
reinforcement.  The  rectangular  planform  is  575  sq.ft 
with  a  wingspan  of  40  ft.  The  2D  airfoil  shape  is  a 
customized  proprietary  design  developed  by  FXC 
Corporation.  The  system  payload  can  be  varied  from 


400  to  600  lbs  without  exceeding  the  performance 
range  of  the  ram-air  canopy.  U  nder  the  current  test 
configuration,  the  payload  is  a  wooden  rectangular  box, 
rigged  to  500  lbs. 

The  guidance  and  control  unit  is  suspended  between  the 
canopy  and  the  payload.  It  consists  of  an  airborne 
computer  with  embedded  control  logic,  GPS  receiver, 

I M  U  and  compass.  T  wo  control  actuators  provide 
control  inputs  to  the  control  lines  attached  to  the  left 
and  right  trailing  edges  of  the  canopy.  The  control 
surfaces  of  the  canopy  can  be  deflected  symmetrically 
(up  to  18  Inches )  for  pitch  control  and  asymmetrically 
for  lateral  control. 


Parafoil  Model 


The  model  was  derived  using  general  equations  of  fluid 
dynamics.21,27 


The  final  form  of  the  equations  of  motion  of  the  whole 
system  as  a  rigid  parafoil-payload  body17,18  was  written 
in  matrix  form 


V*T  =  A  1 
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where  the  following  vector  notation  was  used 
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The  matrices  A,  E  and  U  were  defined  as  follows 
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I  n  these  equations  m  is  the  total  mass  of  the  system, 

/ y  are  the  moments  of  inertia,  cry  are  the  apparent 
mass  and  inertia  terms,  l/,  are  the  velocity 
components,  to,  -  the  angular  velocities,  zs  and  xs  are 
the  vertical  and  longitudinal  moment  arms  from 
particular  component  from  the  origin  of  the  body- 
fixed  coordinate  frame. 

Implementation  was  carried  out  in  the  M  athworks 
M  ATLAB/Simulink  environment  similarly  to  how  it 
was  done  for  the  circular  parachute  in  a  previous 
study.20,21  The  general  view  of  the  Simulink  model  is 
presented  in  Figure  2. 


Aerodynamics  M  odel 

The  only  available  data  on  the  aerodynamic 
coefficients  was  the  wind  tunnel  data  from  Refs.  4 
and  7  and  computational  fluid  dynamics  data 
obtained  forthe  steady  state  turn  in  Ref.  28.  However 
it  was  impossible  to  implement  this  data  directly. 
While  Refs.  7  and  28  contain  the  data  just  for  one 
single  point  in  terms  of  dependence  on  the  angle  of 
attack,  the  data  in  Ref.  4  despite  of  the  richness  of 
different  dependences  for  three  aspect  ratios  is 
somewhat  incomplete  and  contradictory.  Figures  3-7 
give  an  example  of  the  data  on  some  lateral  and 
longitudinal  coefficients  after  some  necessary 
corrections,  it  is  seen,  that  the  data  is  quite  non¬ 
linear,  non-consi stent  and  obviously  use  the  different 
scale  for  the  angle  of  attack.  After  thorough  analysis 
of  the  data  it  was  assumed  that  a  setup  for  the  wind 
tunnel  experiment  included  about  10e  rigging  angle 
so  that  the  angle  of  attack  axis  for  the  longitudinal 
channel  (F  igures  6  and  7)  should  be  shifted  left  by 
this  angle. 

Instead  of  modeling  the  aerodynamic  coefficients 
along  the  full  range  of  angles  of  attack  ([-10  e,80  e])  it 
was  suggested  to  replace  them  with  the  linear 
dependences  within  the  operable  range  of  angles  of 
attack  ([0  e,20  9])  as  shown  by  appropriate  formulas 
on  Figures  3-7.  Not  only  this  provides  a  reasonable 
behavior  and  gives  an  acceptable  accuracy  but  it  also 
enables  the  use  simple  PI  D  technique  where  only  a 
few  parameters  for  each  dependence  have  to  be 
optimized. 

After  careful  analysis  and  preliminary  simulation 
runs  the  applied  aerodynamics  model  was  developed 
as  the  database  for  all  major  aerodynamic  coefficients 
expressed  as  functions  of  several  input  parameters. 
The  longitudinal  and  lateral-directional  aerodynamic 
coefficients  depend  on  the  angle  of  attack  (a), 
sideslip  angle  (0),  symmetrical  flaps  deflection 
(Sf  =  min  (Sleft,  Sright)  e  [0;l]),  and 
differential  flaps  deflection  ( Sa  =  Sleft  -  8right)  as 
shown  in  the  following  equations. 


Fig.  2  Simulink  model 


CD  =  0. 14  +  (0.25  +  0.2Sf)  Cl  (8) 

Cy  =  (-0.005  -  0.0001a) 0  +  (-0.007  +  0 .0012a)  Sa  (9) 
CL  =  0 . 0375  («  +  10°)  +  0 . 1  (2 Sf  +  |da|)  (10) 
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Cn  =  (0.007  -  0.0003a)  J3  +  (0.019  -  0.0008a)  5a  +  (0.02  3 p  -  0.0936 r)  (13) 


whereCt,  CD,  and  CY  are  the  lift,  drag  and  side-force 
coefficients;  C/,  Cm,  and  C„  are  the  rolling  moment, 
pitch  moment,  and  yawing  moment  coefficients;  p,  q, 
and  r  are  the  rolling,  pitch,  and  yawing  moments;  b 
and  c  are  the  inflated  span  and  chord;  V  is  the  total 
velocity;  m1  and  m2  are  the  payload  and  canopy 


masses,  respectively.  The  rigging  angle  of  8e  and 
ac  =0  =  -10°  were  applied. 

The  derivation  of  the  formula  for  the  coefficient  Cm 
was  based  on  the  analysis  presented  in  Ref.  19. 


Alpha  (deg) 

Fig.  4  Modeling  of  the  coefficient  CMa 
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Fig.  5  Modeling  of  the  coefficient  CnSa 


Alpha  (deg) 

Fig.  6  M  odeling  of  the  coefficient  CL 


Alpha  (deg) 

Fig.  7  M  odeling  of  the  coefficient  C D 
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Apparent  Mass  Model 

In  general  the  apparent  mass  tensor  A  (5)  contains  36 
elements  that  maybe  unique  in  a  real  fluid.  For  ideal 
fluid,  however,  the  tensor  /\  has  a  symmetrical  form, 
leaving  a  maximum  of  21  distinct  terms.  In  the  case 
of  a  body  with  one  plane  of  symmetry,  tensor  A  can 
be  further  reduced  to  12  unique  components.21,29,30 

During  descent  of  a  fully  deployed  symmetric 
parafoil  canopy,  its  added  mass  terms  are  functions  of 
air  density  and  are  defined  to  be  the  algebraic 
functions  of  cylindrical  volume  and  parafoil 
geometry  multiplied  by  empirical  scaling 
coefficients.24  Off-diagonal  terms  in  (5)  are 
introduced  as  corresponding  values  of  static  and 
inertial  moments  multiplied  by  correction 
coefficients.  These  unknown  so  far  coefficients, 
which  by  the  way  depend  on  the  state  of  control 
surfaces,  are  left  in  (5)  for  the  further  PI  D . 


Preliminary  Simulation  Results 

The  model  response  on  the  control  inputs  is  discussed 
in  this  section. 

Figure  8  indicates  the  longitudinal  response  to  a  pulse 
input  of  8f  =  o .  5  after  the  parafoil  was  trimmed  in 
flight  at  an  altitude  of  3000  m.  With  flap  deflection, 
the  angle  of  attack  increases,  pitch  angle  increases, 
and  flight-path  angle  is  steeper,  as  is  expected.  A 
short-period  type  of  dynamic  response  is  also  noted, 
which  is  heavily  damped. 

Figure  9  indicates  the  trajectory  due  to  the 
asymmetric  input  of  an  aileron.  Figure  10  shows  roll- 
rate  and  yaw-rate  responses  due  to  the  control  input. 
The  steady-state  bank  angle  is  about  4  2  and  the 
steady  yaw  rate  is  2.8  2/s. 


Fig.  8  Longitudinal  response  to  flap  deflection 
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T rajectory  Data  Source 

The  following  airborne  suite  of  sensors  has  been  used 
to  obtain  the  objective  information  during  this  study: 

•  IMU  block  provided  the  following  data  sampled 
at  100  Hz:  Roll  Rate,  Pitch  Rate,  Y  aw  Rate,  X 
Acceleration,  Y  Acceleration,  Z  Acceleration; 

•  Compass  provided  attitude  data  measured  in  the 
body  frame  and  sampled  at  4  Hz; 

•  GPS  unit  provided  the  following  data  sampled  at 
10  Hz:  Latitude,  Longitude,  Altitude,  Ground 


Speed,  Track  Angle,  Velocity  East,  Velocity 
North,  Velocity  Up. 

The  W  indPack  dropsonde,  dropped  at  the  same  time 
as  the  descending  system,  provided  the  real-time 
wind  profile  update  with  the  frequency  of  10  Hz.  A 
sample  of  vertical  profiles  for  both  parafoil  and 
WindPack  dropsonde  is  shown  on  Figure  11. 

The  data  collected  during  the  only  (so  far)  flight  test 
was  analyzed  and  used  as  a  truth  data  for  the  PID . 


Fig.  11  Example  of  vertical  profiles  versus  time 


Parameter  Identification  Technique 

For  the  case  of  the  rigid  wing  one  of  the  well- 
established  system  ID  tools  may  be  used  to  identify 
its  aerodynamic  characteristics.  In  the  parafoil  case 
however  this  approach  is  impossible.  Moreover,  the 
estimation  of  the  apparent  mass  terms  even  for  the 
case  of  the  fixed-geometry  parafoil  is  non- 
i ntuiti ve.21,28,29  T he  apparent  mass  has  the  following 
important  property:  it  changes  with  orientation  of  the 
associated  motion.  Therefore,  appropriate  choice  of 


truth  data  can  obviously  identify  the  apparent  mass 
terms  that  dominate  specific  motion.  A  non-standard 
nonlinear  system  PID  technique  including 
employment  of  the  multicriteria  optimization31  to  find 
reasonable  ranges  of  the  optimization  parameters  and 
zero-order  FI  ooke-J  eeves  method32  to  actual  ly  def i  ne 
the  best  solutions  for  optimization  parameters  was 
applied  to  tune  the  initial  aerodynamic  dependences 
and  apparent  mass  terms  as  well.  Although  the  PID 
still  continues,  Figures  12-14  show  some  preliminary 
results  obtained  so  far. 


7 


Horizontal  and  Vertical  Projections 
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Fig.  12  Comparison  of  the  flight  test  data  with  the  original  (not-tuned)  model 


Horizontal  and  Vertical  Projections 


2000 


£  1500 

CD 

T3 

13 

§j  1000 
< 


500 


Fig.  13  Comparison  of  the  flight  test  data  with  the  tuned  model 
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While  the  original  model  matches  the  integral 
behavior  of  the  system  in  terms  of  glide  ratio  (2.5) 
and  vertical  velocity  (3.7  m/s)  fairly  well  some 
correction  was  needed  to  match  the  turn  rate  6  °/s 
demonstrated  during  the  flight  tests. 

Figure  12  demonstrates  the  discrepancy  of  the 
original  model  in  the  lateral  channel.  Figure  13  shows 


how  the  changing  of  the  only  coefficient  from 
0.019  to  0.032  i  mproves  the  model .  F  urther 
improvement  involving  more  coefficients  is  also 
possible  (Figure  14).  Flowever  in  this  case  the 
developed  PID  technique  involves  finding  and  using 
only  Pareto-optimal31  solutions,  which  do  not  make 
other  criteria  (vertical  channel  matching,  eigen 
frequencies,  etc.)  worse. 
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Horizontal  and  Vertical  Projections 


Fig.  14  Exampleof  the  further  tuning 


Conclusions 

The  complete  6D0F  model  of  the  parafoil  system  has 
been  developed  and  realized  within  the  M  athWorks 
M  ATLAB/Simiulink  environment.  The  model 
demonstrates  an  adequate  response  for  the  control 
inputs  and  matches  the  result  of  the  flight  tests  well. 
A  nonlinear  system  identification  technique  is 
currently  being  applied  to  tune  the  model  and  correct 
the  apparent  mass  terms.  The  authors  intend  to 
present  more  PI  D  results  upon  the  completion  of  the 
tests. 
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